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Long  Term  Goals 

•  Develop  a  coupled-domain  formulation  that  will  lead  to  accurate  predictions  of  the  dynamics 
and  structural  integrity  of  F-IX  and  F-22  aircrafts  in  action  under  dynamic  and  extreme  loading 
conditions. 

•  Implement  new  algorithms  for  high  Reynolds  number  flows  in  complex  geometries  as  well  as 
methods  for  the  nonlinear  acroelastic  motion. 

•  Address  issues  of  heterogeneous  discretizations  and  corresponding  software  integration  ap¬ 
proaches  within  the  parallel  environment  context. 

Background  and  Specific  Objectives 

Modeling  of  the  buffeting  and  flutter  regimes  require  both  an  accurate  approximation  of  the 
dynamic  loading  produced  by  the  fluid  as  well  as  a  realistic  model  of  the  structure.  Most  of  the 
previous  works  involve  simplifying  assumptions  either  in  the  fluid  or  the  structural  solver.  On  the 
part  of  the  fluid  solvers,  assumptions  like  using  the  inviscid  equations  or  ad  hoc  turbulence  models 
are  used  to  simplify  the  fluid  computation.  On  the  part  of  the  structure,  thin  solid  structures  are 
modeled  with  combinations  of  beam,  plate  and  shell  elements,  each  of  which  are  formulated  upon 
certain  assumptions. 

Under  this  AFOSR  grant,  we  have  developed  a  parallel  spectral/hp  element  fluid  solver  that: 

•  enables  an  accurate  and  fast  resolution  of  the  flow  field,  and 

•  is  linked  to  an  hp-FEM  structural  solver,  which  in  turn  enables  a  realistic  representation  of 
thin  solid  structures,  e.g.  wings. 
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Specifically,  the  flow  solver  is  AfsKTocr,  which  we  have  developed  and  modified  for  LES  type 
computations  of  compressible  flows.  The  structural  solver  is  StressCheck,  which  has  also  been 
initially  developed  under  AFOSR  support  by  Prof.  B.  Szabo  and  subsequently  by  Engineering 
Software  Research  &  Development,  Inc.. 

Platform  Independence  Strategy  for  Fluid-Structure  Coupling 

Two-way  coupling  of  MsKTotr  and  StressCheck  requires  a  software  solution  for  extracting 
traction  information  from  the  fluid  solver  and  providing  it  to  StressCheck.  At  this  stage,  a 
non-linear  elastic  analysis  is  performed,  and  the  deformations/ velocities  on  the  structure  faces  are 
passed  back  to  A/ewT" otr  for  updating  the  position  of  the  structure  within  the  fluid.  This  can  be 
achieved  by  software  integration  in  such  a  way  that  both  codes  could  interact  dynamically. 

Because  MskTolt  runs  under  the  Unix  operating  system,  and  StressCheck  runs  under  the 
Microsoft  Windows  environment  a  software  solution  had  to  be  devised  to  overcome  this  obstacle. 
The  software  solution  we  have  chosen  is  to  connect  the  two  software  entities  via  sockets.  The 
procedure  goes  as  follows.  Using  StressCheck’s  Component  Object  Model  (COM)  interface,  a 
Visual  Basic  socket  server  application  was  developed  which  listens  on  a  specified  port  on  the  PC 
on  which  the  application  is  running.  Once  the  StressCheck  server  is  initiated,  the  client  (in  the 
sense  of  socket  communication)  fluid  solver  is  executed  on  a  parallel  Unix  platform,  and  attempts 
connection  to  the  server  PC  on  the  pre-specified  port.  Once  a  socket  connection  is  initiated,  data 
may  be  transferred  between  the  two  applications.  A  schematic  diagram  of  the  computer  code 
coupling  is  provided  in  Figure  1. 

Tests  indicate  that  transferring  via  sockets  2.5MBytes  (containing  the  pressure  field  on  a  wing 
as  obtained  by  AfsKT aron  a  Unix  platform  and  transferred  to  StressCheck  on  a  PC  platform) 
lasts  less  than  1  second. 

Prototype  Problems 

Two  model  problems  are  considered  for  the  verification  and  validation  of  the  fluid-structure 
coupling:  the  AGARD  wing  and  the  square  flat  plate  subjected  to  flow. 

The  AGARD  wing: 

We  considered  the  semi-span  solid  wall-mounted  wing  with  a  quarter-chord  sweep  of  45°  and 
a  NACA  65A004  aerofoil  section.  The  semi-span  of  the  wing  is  s  =  2.5ft,  the  root  chord  is 
tst  =  1.833/<,  and  the  tip  chord  is  tbt  =  1.208/t.  In  Yates  et.  al  [1]  two  wings  of  having  the 
same  dimensions  were  wind-tunnel  tested.  A  hp-finite  element  parametric  model  {s,  tbs,  tbt  can  be 
changed  and  the  overall  finite  element  mesh  changes  accordingly)  consisting  of  48  hexahedral  p- 
elements  and  12  pentahedral  p-elements  was  constructed  using  the  software  product  StressCheck, 
as  shown  in  Figure  2.  The  structmal  finite  element  mesh  is  refined  in  the  neighborhood  of  the  lead¬ 
ing  edge  to  better  represent  the  pressure  sharp  gradients,  and  at  the  root  to  control  the  numerical 
error  due  to  the  singularities  associated  with  the  clamped  boundary  conditions. 

The  total  mass  of  the  wing  as  reported  in  Yates  et/  al[l]  is  0.145025  ±0.001555  slugs,  and  its 
volume  as  computed  by  the  FE  model  is  0.159658  ft^.  These  two  values  are  used  to  compute  the 
material’s  density:  ps  =  0.908344  ±0.0097  slug/ft^. 


2 


Figure  1:  ^e/cT'ar-StressCheck  coupling  strategy  via  sockets. 


As  the  wing  is  made  of  mahogany  wood,  it  is  an  orthotropic  material  with  mechanical  properties 
which  are  uncertain.  However,  as  claimed  by  the  author  who  performed  the  tests  in  Yates[2],  the 
material  is  homogeneous.  Its  mechanical  properties  can  be  found  in  two  references[4]  and  [5],  which 
are  summaries  of  data  taken  from  Anon[6]  (cited  by  Yates[2]).  The  fiber  orientation  of  the  wood 
is  taken  along  the  wing  span,  which  is  inclined  45  degrees  with  respect  to  y  axis.  Based  on  the 
references  above,  we  summarize  in  Table  1  the  diflFerent  mechanical  properties  of  the  mahogany 
wood.  Eii,E22  and  E33  denote  the  Young  modulus  in  the  longitudinal  (fiber),  tangential  and 
radial  direction,  respectively. 

Validating  the  CSM  AGARD  wing  model 

To  check  the  validity  of  the  results  obtained  by  the  two-way  coupling,  we  first  validated  the 
computational  structural  model  (CSM).  This  was  accomplished  by  computing  the  eigen-frequencies 
as  obtained  by  the  numerical  simulation  and  comparing  them  with  the  experimental  observations. 
Hence  we  assume  that  both  numerical  and  modeling  errors  of  the  structural  simulation  are  negli- 
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Figure  2:  Computation  fluids  mesh  and  top  view  of  the  computational  mechanics  mesh  (top)  and 
dimensions  for  the  AGARD  problem  (lower). 


gible. 

As  observed,  the  Mahogany.l  and  Mahogany _4  are  very  similar  in  terms  of  their  mechanical 
properties.  We  performed  hp-FEM  eigen-analysis  using  the  three  different  material  properties, 
increasing  the  p-level  from  1  to  6  (at  p  =  6  there  are  8460  degrees  of  freedom),  with  a  numerical 
error  in  the  flrst  flve  frequencies  computed  by  a  posteriori  error  estimators  being  less  than  0.005%. 
A  t5q)ical  convergence  curve  of  the  5*^  frequency  is  shown  in  Figme  3.  This  graph  shows  that  the 
numerical  errors  are  negligible. 

In  Table  2  we  summarize  the  first  five  eigen-frequencies  of  the  wing  obtained  by  om:  FE  model, 
the  ones  computed  in  Yates[2]  and  the  experimental  observations  in  the  last  two  colunms.  In 
Yates[2]  a  2-D  approximation  has  been  applied  by  assuming  plate  elements  and  material  properties 
which  have  been  altered  (by  more  than  10%)  in  order  to  better  “approximate”  the  test  results.  Note 
that  using  the  hp-FEM  structural  code,  no  such  assumptions  are  necessary. 

One  may  notice  the  good  agreement  between  the  hp-FEM  results  using  material  properties 
Mahogany.l  and  the  test  results.  This  assures  that  our  model  is  free  of  both  numerical  and  modeling 
errors  (a  proper  3-D  anisotropic  representation  of  the  wing,  and  material  properties  Mahogany_l 
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Table  1:  Mechanical  properties  of  the  laminated  mahogany 


Symbol 

Mahogany -1 
Ref.  [5] 

Mahogany_4  “ 
Ref.  [4] 

Ell  10“  Ibf/fP 

192.96 

192.96 

E22  10®  Ibf/fP 

12.63 

12.35 

E33  10®  Ibf/fP 

21.11 

20.65 

Gi2  10®  Ibf/ft^ 

13.02 

16.59 

Gi3  10®  Ibf/fP 

16.59 

12.74 

G23  10®  Ibf/fP 

5.53 

5.40 

vn 

0.034 

0.533 

V13 

0.033 

0.314 

l'23 

0.326 

0.326 

Honduras  Mahogany  (Swietenia  macrophylla) 


in  “Green”  moisture  content. 


Table  2: 

First  five  modal  frequencies  I 

(Hz) 

# 

Mah_l 

Mah_4 

Comp.  [2] 

Exp.  [1] 

Exp.  [1] 

1 

14.37 

14.53 

14.12 

14.60 

14.10 

2 

47.99 

52.43 

50.91 

47.70 

50.70 

3 

69.94 

71.04 

68.94 

67.30 

69.30 

4 

119.01 

127.76 

122.25 

117.00 

127.10 

5 

169.51 

171.69 

160.52 

used  based  on  literature  reference  without  alteration). 

The  necessity  of  geometrically  non-linear  analysis  of  wings  is  illustrated  by  considering  the 
AGARX)  wing  at  a  20"  AOA.  In  Figure  4  we  present  the  principal  stresses  on  the  deformed  shape 
as  the  wing  is  loaded  by  the  pressure  field  obtained  by  AfenT" ar.  We  note  the  large  difference  in 
the  deflected  geometry  and  the  stress  level  between  the  linear  and  non-linear  analyses. 

Linear  vis  Non-linear  Structural  Analysis  and  Pressure-  Velocity  Coupling 

Because  the  finite  element  mesh  for  the  computational  fluids  simulation  is  completely  different 
compared  to  the  mesh  suitable  for  the  solid  mechanics  simulation,  two  different  meshes  are  used. 
For  example,  the  AGARD  wing,  modelled  as  a  fully  3-D  sohd  in  StressCheck  has  “needle  type” 
elements  with  aspect  ratios  close  to  1000,  as  opposed  to  the  mesh  used  for  the  fluid  simulation  ( 
Figure  2).  Thus,  the  pressure  field  passed  from  AfeK'Tar  to  StressCheck  and  the  velocities  due 
to  the  deflection  of  the  wing  structure  passed  back  from  StressCheck  to  A/eK'Tar,  have  to  be 
correctly  transferred  as  to  not  destroy  the  accuracy  properties  of  either  the  fluid  or  the  structural 
simulation.  Because  high-order  schemes  and  blending  function  mapping  methods  are  used,  we  can 
specify  these  as  functions  on  the  surfaces  which  are  represented  exactly  (including  the  curvatures) 
in  both  codes.  The  pressure  fleld  is  extracted  on  four  surfaces  on  the  wing  in  AfsK'Tcxr,  the 
upper  and  lower  surfaces,  leading  edge  and  tip  surface.  For  the  AGARD  wing  about  10,000  points 
and  their  pressure  values  are  passed  through  sockets  to  the  VB  COM  interface.  Each  of  the  four 
surfaces  is  projected  onto  a  plane,  and  a  pressure  functional  representation  is  constructed  using 
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Emx  EsL:  SoUion  ••  SOLjbObiM.  runs  tti  to  tiS  (^L  ^  Elements) 


Estimated  LM  1.19003423e^[cycles/sec) 
Estineted  Limit  7.47720598e+002  (ladartt/sec) 
NOTE:  If  time  not  seconds,  leedfcvdes/limel  ot  fjadians/timel 


Figure  3:  Typical  convergence  rate  in  the  computation  of  eigen-frequencies.  The  5*^  frequency 
versus  DOF  for  Mahogany _1  material. 


least-squares  of  iV-th  order.  For  example,  the  upper  and  lower  wing  surfaces  are  projected  onto 
the  X  —  y  plane,  and  the  pressure  is  represented  by: 

N 

i=0 

=  Uo  +  0,iX  +  02?/  +  0-ZX^  +  0-AXy  +  05?/^  H - 

The  number  of  terms  in  the  pressure  approximation  series  is  determined  by  an  adaptive  pro¬ 
cedure,  namely,  the  functional  representation  of  the  the  pressure  over  each  surface  is  computed 
for  increasingly  larger  iVs,  and  the  total  force  is  extracted.  At  a  given  N  the  total  force  remains 
unchanged,  so  this  is  the  value  used  in  (1)  for  the  representation  of  the  pressure.  Preliminary  tests 
have  shown  that  good  results  are  obtained  for  the  problems  of  interest  for  N  >  10,  and  the  total 
force  computed  through  the  use  of  (1)  in  StressCheck  is  almost  identical  to  the  total  force  on  the 
wing  as  computed  at  the  10,000  points  in  A/sKTar. 

To  demonstrate  the  need  for  a  structural  geometric  non-linear  analysis  to  realistically  represent 
the  wing  in  flight,  the  AGARD  wing  has  been  analyzed  as  the  angle  of  attack  is  increased  up  to 
55°.  The  data  has  been  transferred  to  StressCheck  together  with  the  data  poo  =  1/2  Atm  and 
Pinfty  of  stajidard  air.  After  the  fluid  solver  passes  the  pressure  information  to  the  StressCheck 
server  through  the  socket  interface,  the  later  computes  the  displacements  and  velocities  due  to  the 
newly  updated  loading.  In  Figure  5  the  structural  response  for  a  linear  and  geometrically  non-linear 
analysis  was  performed  showing  a  big  difference  between  the  two  cases.  It  is  noticed  that  the  more 
realistic  non-linear  analysis  provides  displacements  of  a  much  larger  magnitude  (especially  at  the 
large  angle  of  attack). 
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Figure  4:  The  AGARD  wing  at  angle  of  attack  of  20®,  linear  (left)  vis  non-linear  (right)  responses. 
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AO  A  =  30° 


AO  A  =  55° 


Figure  5:  The  principal  stresses  in  the  AGARD  wing  at  two  angles  of  attack:  linear  (left)  vs 
geometrical  non-linear  (right)  response. 


The  flexible  thin  plate  in  flow 

A  second  benchmark  problem  is  considered  which  is  both  simple  geometrically,  and  does  not 
require  large  computational  resources.  Consider  a  plate  made  of  isotropic  material  of  dimensions 
1  X  1  X  h  where  h  «  1  is  considered,  clamped  at  its  lateral  boundaries.  In  the  computational 
fluids  mesh  the  plate  is  assumed  to  be  part  of  the  x  —  y  plane  where  a  flow  is  present  for  >  0  and 
a  constant  pressure  po  is  applied  on  the  face  z  =  —h.  This  example  problem  is  described  in  detail 
in  Gordnier  and  Visbal  [3],  and  was  used  as  a  benchmark  problem.  The  CSM  finite  element  model 
consists  of  128  hexahedral  p-elements,  with  2  elements  in  the  thickness  direction  and  64  elements 
in  the  x  —  y  plane  (with  needle  elements  needed  to  accurately  capture  the  solid  boundary  layers). 
The  CSM  and  CFD  meshes  used  for  this  test  problem  are  shown  in  Figure  6. 


The  von-Karman  Plate  Model 


NeKxar: 

Re=10,000 
M  =  0.5 
At  =  0.001 

p=  2"“  (up)  3"^  (down)  order  per 
direction 

720  hexahedra  elements 
interfacing: 

MR:  1000,  RR:  10000 
von-Karman: 

a=b=1,  h=0.01  (used  h=0  since  no 

difference) 

v=0.3 

c=  1.25 


N=13 

At  determined  by  scaling 


Figure  6:  CFD  mesh  and  streamlines  for  the  square  plate  problem. 


We  developed  numerical  algorithms  for  treating  the  von-Karman  plate  model  which  is  a  nonlin¬ 
ear,  dynamic,  partial  diflferential  system  over  rectangular  domains.  It  is  solved  by  the  Chebyshev- 
collocation  method  in  space  and  the  implicit  Newmark-/3  time  marching  scheme  in  time.  In  the 
Newmark-/?  scheme,  a  non-linear  fixed  point  iteration  algorithm  was  employed. 
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We  monitored  both  temporal  and  spatial  discretization  errors  based  on  derived  analytical  so¬ 
lutions,  demonstrating  highly  accurate  approximations.  We  also  quantified  the  infiuence  of  a  com¬ 
mon  modeling  assumption  which  neglects  the  in-plane  inertia  terms  in  the  full  von-Karman  system, 
demonstrating  that  it  is  justified.  A  comparison  of  our  steady-state  von-Karman  solutions  to  pre¬ 
vious  results  in  the  literature  and  to  a  three-dimensional  high-order  finite  element  analysis  was 
performed,  showing  an  excellent  agreement.  Other  modeling  assumptions  such  as  neglecting  in¬ 
plane  quadratic  terms  in  the  strain  expressions  were  also  addressed  -  see  details  in  the  publications 
list. 

The  two-way  coupling  algorithm  has  been  implemented  to  simulate  the  von-Karman  plate  in 
flow,  where  the  developed  von-Kdrman  code  has  been  coupled  with  Afen'T arvia  sockets.  Figure 
7  shows  the  PSI  analysis,  monitoring  the  mid-plane  displacement,  and  convergence  in  terms  of 
time-step  parameter. 


von-Karman  -  Nektar  coupUng 
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f  1  «-r< 


«  !  ; 


tc 


1  imr 


ft,  M 


Figure  7;  AfsKTar-  von-Karman  coupling  for  simulation  of  a  square  plate  undegrgoing  non-linear 
deflection  due  to  the  flow  field. 
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It  is  also  demonstrated  that  to  ensure  control  of  the  discretization  error  and  that  the  flow  field 
has  been  resolved,  the  p— level  over  the  elements  has  been  increased  until  convergence  has  been 
obtained.  This  influence  is  shown  in  Figure  8,  which  demonstrate  that  the  mid-plane  transverse 
displacement  requires  the  p-level  to  increase  up  to  p  =  5  to  obtain  converged  results. 


Figure  8:  A/"e«T"ar-  von-Karman  coupling:  Convergence  in  the  response  as  p-level  is  increased. 


Remaining  Open  Questions 

At  the  moment  StressCheckdoes  not  enable  solutions  of  time  dependent  problems.  A  simplistic 
approach  has  been  incorporated  to  account  for  inertial  terms  (due  to  acceleration)  which  performs 
well  only  for  small  acceleration,  and  diverge  for  large  acceleration.  For  a  successful  FSI  scheme, 
time  dependent  capabilities  for  finite  deformations  have  to  be  developed  and  implemented  into 
StressCheck.  These  activities  had  started  and  are  based  on  F.  Armero  research  which  was  also 
supported  by  the  AFOSR.  This  newly  created  fluid-structure  hp-based  code  will  also  be  used  in 
the  future  for  addressing  questions  related  to  survivability  of  aircraft  when  parts  of  the  structure 
axe  hit  or  torn. 
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